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Abstract 

A new framework for analysing the gravitational fields in a stationary, axisymmetric 
configuration is introduced. The method is used to construct a complete set of field 
equations for the vacuum region outside a rotating source. These equations are under¬ 
determined. Restricting the Weyl tensor to type D produces a set of equations which 
can be solved, and a range of new techniques are introduced to simplify the problem. 
Imposing the further condition that the solution is asymptotically flat yields the Kerr 
solution uniquely. The implications of this result for the no-hair theorem are discussed. 
The techniques developed here have many other applications, which are described in 
the conclusions. 


1 Introduction 

The field equations of general relativity are remarkable both for their complexity and for 
the incredible variety of techniques that can be applied in generating solutions. These 
range from brute force computation through to advanced techniques based on symmetries 
and invariants B0- In this paper we introduce a new method for tackling the gravi¬ 
tational field equations for stationary, axisymmetric systems. The main applications of 
this scheme are to finding the gravitational fields inside and outside rotating sources. The 
method has strong similarities with the Newman-Penrose (NP) formalism Q, and also 
with coframe methods for analysing differential equations. The essential difference with 
the NP formalism is that our approach is based on a real orthonormal tetrad, as opposed 
to a complex null tetrad. The only complex structure we work with arises naturally in the 
structure of the bivector fields generated by the tetrad vectors. 

The main advance in the work presented here is that the technique can be pushed 
right through to an explicit solution. This is possible with the introduction of a number 
of new methods. The closest precursor to these techniques is the work of Held ||], though 
we believe our current methods offer many improvements. In this paper we concentrate 
on the vacuum equations. We extract a complete set of vacuum equations, and find that 
these are under-determined for a general, axisymmetric source. To construct a unique 
solution we then impose the additional condition that the Riemann tensor is type D, and 
the solution is asymptotically flat. We are then able to show constructively that the Kerr 
is the unique solution for this case. This result does not rely on the existence of a horizon. 
The Carter-Robinson uniqueness theorem can therefore be interpreted as stating that the 
formation of a horizon is related to a restriction on the algebraic form of the Weyl tensor. 
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The Kerr solution occupies a unique place amongst the known exact solutions to the 
Einstein field equations. Many derivations of the solution have been found, a number 
of which are summarised in Chandrasekhar’s classic text Few authors (apart from 
Chandrasekhar) have claimed that the derivation of the Kerr solution is simple, and many 
of the derivations found are highly opaque. Probably the most mysterious derivation of all 
is that of the complex coordinate transformation trick [j|, In this derivation one starts 
with the Schwarzschild metric in advanced Eddington-Finkelstein coordinates, expressed 
in terms of a complex null tetrad. A complex coordinate transformation is then applied 
to yield a new metric, which turns out to be that of the Kerr solution. This is a trick 
because there is no reason to expect the complex coordinate transformation to generate a 
new vacuum solution. The justification for this is only revealed at the end of a detailed 
study of vacuum metrics of Kerr-Schild type, and is quite obscure |]] . 

The derivation presented here proceeds in a very different way to the usual ‘metric’ 
route adopted in general relativity (GR). In a typical GR derivation one starts with a 
metric of suitable symmetry containing a set of arbitrary scalar functions. The Einstein 
equations are then written as a set of coupled second-order differential equations in the 
metric functions. These equations are notoriously difficult to analyse, for many different 
reasons. One of the difficulties is that the freedom to perform coordinate transformations 
must be removed by restricting the form of the metric (‘gauge-fixing’). If this is not done 
correctly, one has more unknowns than there are equations for. The problem is that it is 
often not clear ab initio how best to perform this gauge fixing. Ideally one would like the 
chosen coordinates to have some simple physical interpretation, but how to achieve this 
does not usually emerge until later in the solution process. 

Our method avoids this problem by developing a first-order set of equations analogous 
to those obtained in the NP formalism. Gauge-fixing of the Lorentz group is performed 
at the level of the Riemann tensor, which makes it easier to ensure that the gauge choices 
are physically sensible and mathematically convenient. The resulting equations relate 
abstract derivatives of the terms in the spin connection to quadratic combinations of the 
same quantities. This provides a very clear way of expressing and analysing the non- 
linearities in the theory. Once the Riemann tensor has been found, some natural physical 
scalar fields start to emerge. Invariance under diffeomorphisms is then employed to give 
these physical fields simple expressions in terms of our chosen coordinates. This process 
lifts the status of the coordinates from arbitrary spacetime functions to local physically- 
measurable fields. We then introduce a series of new techniques which solve directly for 
the spin coefficients without computing any of the metric coefficients. The latter are then 
found by straightforward integration. This is a considerable advance on previous first-order 
formulations of GR, and opens up a number of new solution strategies. We do not yet 
possess a general method for solving all problems in our new formalism, but there are good 
reasons to believe that some very general techniques do exist. The derivation of the Kerr 
solution presented here provides a number of clues to how such general methods might be 
found. 

One feature of our method, which it shares with the NP formalism, is that the early 
stages involve simple, repetitive algorithms for organising terms. This work is well suited 
to symbolic algebra packages such as Maple. The derivation of the equations is further 
simplified by employing the language of Clifford algebra i I 0- 0- This is employed 
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as a set of algebraic rules for manipulating vectors, rather than via a concrete matrix 
representation (such as the Dirac matrices). In this sense we are adopting the language 
of spacetime algebra [12, [H|. This approach to Clifford algebra is also well suited to 


implementation in symbolic algebra packages. The application of spacetime algebra clearly 
demonstrates the origin of a natural complex structure generated by the spacetime directed 
volume element, or pseudoscalar. This element squares to minus one, so unites the complex 
structure of the 2-spinor approach with duality transformations in tensor language. This 
is a considerable unification, and provides a number of insights into the complex structure 
underlying the Kerr solution. This also sheds some light on the origin of the mysterious 
complex coordinate transformation trick discovered by Newman and Janis [5|]. 

This paper starts by introducing the variables we use to encode a stationary axisyrn- 
metric gravitating system. We then introduce a suitably general ansatz to describe axially- 
symmetric fields. The link between the gauge fields and the metric and Christoffel connec¬ 
tion is explained. The complete set of vacuum equations is obtained, and we then restrict 
to looking for vacuum solutions of type D. Imposing this restriction causes the equations to 
simplify in a truly remarkable manner to a set of 8 coupled equations. These equations are 
then simplified by identifying various integrating factors. This enables us to make a range 
of deductions about the solution without introducing a definite coordinatisation. Finally, 
a concrete expression of the solution is produced, which reproduces the line element of the 
Kerr solution in Boyer-Lindquist coordinates. We end with a discussion of the insights the 
present formalism brings to the no hair theorem for black holes. This states that the the 
Kerr solution is the only possible vacuum axisymmetric solution outside a horizon, which 
is now interpreted in terms of a restriction on the algebraic type of the Weyl tensor at a 
horizon. 

Summation convention and natural units (c = G = 1) are employed throughout. Greek 
letters are employed for coordinate indices and Latin for tetrad components. We work in 
a space with signature (1, —1, —1, —1). 


2 Field equations for axisymmetric systems 


The equations derived here follow the route developed in an earlier paper [ 101, though 
the present derivation is streamlined and much of the less conventional terminology has 
been removed. A small price for this is that some results in the following derivation are 
stated without proof. Further details can be found in |Io| , 131. Our starting point is the 


Clifford algebra of Minkowski spacetime. The (constant) generators of this are denoted by 
(7o • • ■ 73 } and satisfy 


lilj + 7j7* = i,j = 0...3, (1) 

where our convention is that r/ t j = diag(l, —1, —1, —1). We suppress any mention of the 
identity matrix, which is superfluous for all calculations. Throughout we use Latin indices 
for tetrad frames and components, and Greek for coordinate indices. At various points we 
adopt the language of spacetime algebra and refer to the Clifford generators as ‘vectors’. 
The symmetric and antisymmetric parts of the Clifford product of two vectors define 
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the inner and outer products, and are denoted with a dot and a wedge respectively: 


li'lj = + 7j7i) = Vij 

= ( 2 ) 


A full 

(real) basis for 

the Clifford algebra is provided by 




1 

{7;} 

{7»A7j} {Ty*} 

I 



1 scalar 

4 vectors 

6 bivectors 4 trivectors 

1 pseudoscalar, 

(3) 


grade 0 

grade 1 

grade 2 grade 3 

grade 4 


where 



I = 707x7273. 


(4) 


The 4-vector I is the highest grade element in the algebra and is usually called the pseu¬ 
doscalar, though directed volume element is perhaps more appropriate. The pseudoscalar 
squares to — 1 , 

I 2 = -1, (5) 

and I anticommutes with all odd-grade elements and commutes with even-grade elements. 

We next introduce a set of generators appropriate for the study of axisymnretric fields. 
With (r, 6 , (f>) denoting a standard set of polar coordinates we define 

7 1 = To 

7 r = sin0(cos07i + sin</>72) + cos@73 
70 = cos0(cos07i + shn/> 72 ) — sin $73 

70 = - sin07i + cos072. (6) 

These vectors also form an orthonormal basis for Minkowski spacetime. Henceforth the 

Latin indices i = 0...3 are assumed to run over the set { 7 t, 7 r , 7 e, 70 }- The reciprocal 

frame is given by 

7* = 7 1, Y = -7r, 7^ = -79, 7 0 = -70- (7) 

The (Clifford) product of two orthogonal vectors results in a bivector. There are six of 
these in total, and a useful basis for these is provided by 


Or — 7r7i 


c 

II 

b 


°<t> = 7</>7*- 

(8) 

These all have unit square, 


(T r (J r — (JqCTq — (7'0U'0 — 1 

(9) 

and distinct bivectors anticommute, 


cr (T q — ct q cr , etc. 

(10) 
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The three bivectors also satisfy the identity 

O r <76(7<t> = ItlrlOl# = 70717273 = I- (H) 

The bivector basis is completed by the dual bivectors 

I Or = -lelt 
log = 

lOf = ~7r7f?- ( 12 ) 

The algebra generated by the set {ay, <jq, 07 ,} is isomorphic to the Pauli algebra. This is the 
reason for the notation. It should be remembered, however, that each of the {ay, ag, o^} 
set anticommutes with 7 j. 

The space of bivectors generates the Lie algebra of the group of Lorentz transforma¬ 
tions. Under a proper orthochronous Lorentz transformation a vector a is mapped to the 
vector a 1 according to the rule 

a >—>■ a' = RaR (13) 

where R is an even-grade Clifford element satisfying 

RR = 1. (14) 

The tilde on R denote the reverse operation, which reverses the order of all products in a 
general element. The effect of this on the grade-r element M r is 

M r = (—l) r ( r- 1 )/ 2 M r . (15) 

The object R is referred to as a rotor. Rotors form a group under the Clifford product, 
which gives a spin-1/2 representation of the proper orthochronous Lorentz group. All 
rotors can be written as ± exp(5/2), where B is a bivector. The commutator of a bivector 
and a grade-r object results in a new object of grade-r. For this operation it is useful to 
define the commutator product by 

MxN = \(MN - NM), (16) 

where M and N are general Clifford elements. The space of bivectors is closed under the 
commutator product, and this generates the Lie algebra for the rotor group. 

The key to our approach is to parameterise the gravitational fields via a coframe. We 
introduce the four vector fields 


9 t = 017 * + 

r suit/ 

g r = hY + —7 e 

r 

g e = —j e + C27 r 

V 

9* = —+ 

r suit/ 


(17) 
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where all of the variables (ai ... d 2 ) are scalar functions of r and 6. The indices on the 
{g f > 9 r i g e ) 9^} are to be read as coordinate indices, and we denote the set by { g M }, g = 
0... 3. The reciprocal frame { g to the { g frame is dehned by the equation 

^ + 5^ = 2^. (18) 

The g^ and 7 * vectors dehne a tetrad field in the obvious manner, 

9u,Y +l l 9n = 2V- ( 19 ) 

The metric tensor g^ is derived from the g^ vectors via 

9[i9v + 9v9ii = ^9iiv (20) 


This results in the line element 


, 9 d\ 2 — do 2 r 2 sin 2 # , 9 a\d^r 2 sin 2 # — diOo , ,, ai 2 r 2 sin 2 

ds 2 = =- rr^dt 2 + V , -rTTr^ 2 dt # - 


a 2 


(ai#i - a 2 d 2 ) 2 
ci 2 + c 2 2 r 2 
(ci&i - c 2 b 2 ) 2 


(aidi - o 2 d 2 ) 2 
^ + 6,^ _ 
(C 1 &! - c 2 b 2 ) z 


(aidi - a 2 d 2 ) 2 

61 V + 6 2 2 


■dcj) 2 


(cih - c 2 b- 2 ) 2 


d#" 


( 21 ) 


At first sight this looks unnecessarily complicated. Various nonlinearities are introduced, 
and there are redundant degrees of freedom in the g^ vectors. But our method works 
directly with the g^ and not with the line element. This gives greater control over the 
nonlinearities and the extra (gauge) degrees of freedom can be employed to advantage later 
in the calculation. 

Our main differential operators are the coframe derivatives L,, 


Li = 1%-g^d^ ( 22 ) 

where we employ the abbreviation 

, _ _d_ 

11 dx^ 

The coframe derivatives can be viewed as directional derivatives, with directions deter¬ 
mined by the coframe. Some texts prefer to use a partial-derivative notation for coframe 
derivatives. We have not chosen such a notation because coframe derivatives do not nec¬ 
essarily commute. Written out explicitly the Lj operators are 



L t = aid t + d 2 d<p 
L r = b\d r + c 2 dg 

L g = -{cidg + b 2 d r ) 
r 

L< f> = ^§^ d <t> + a ^ d t)- ( 24 ) 

The coframe derivatives satisfy a set of bracket identities 

[Li,Lj\ = dj k Lk- (25) 


6 










The spin connection is encoded in a set of bivector fields io t . These are related to the Cij k 
coefficients by 

°ij k = l k ■ {Lilj + Wj x 7 j - L^i - Uj x 7 i). (26) 

This equation embodies the content of the first structure equation. 

Rather than solving this equation directly, our method involves introducing a general 
parameterisation for the connection terms. An appropriate parameterisation is defined by 


uj t = -(T + IJ)a r - (S + IK)(Jq - d 2 7172 
i 0 r = —(S + IK)Io + c 2 Ia$ 

UJq = ~(G + 13)1(7$ + —I(7$ 

co$ = -(H + IK)Ia r + {G + IJ)Iag - -4^7172- (27) 

r sind 


This parameterisation introduces a set of 10 scalar functions (G, G, J, J, S, S, K , K , T, 
H). Each of these is a scalar function of r and 9. The reason for the labelling scheme will 
become apparent when the final set of equations is derived. Equation (26) now produces 
the 6 bracket relations 


[Lt, L r ] = —TL t - (K + K)L$ [L r , L e \ = -SL r - GL e 

[L t , Lg\ = —SL t + (J - J)L$ [L r , L$\ = —{K - K)L t - GL$ 

[Lt-,L$\ = 0 [Lg,L$] = (J + J)Lt — HL$. (28) 

The Christoffel connection coefficients are recovered from our various fields through the 

formula 

=9 X -(d l i9v +9 h-Y ViXg u ). (29) 

Clearly this cannot be computed until the g^ vectors are explicitly constructed, which is 
not achieved until the end of the calculation. 


3 Gauge freedom 


There are two distinct types of gauge freedom present in the setup described above. These 
correspond to the two symmetries on which the gauge treatment of gravity is based [|I(], 14]. 
The first is invariance under diffeomorphisms. This can be encoded in various different 
ways, but for our present purposes it is simplest to consider two new scalar coordinates 
r'(r,0) and O'(r,0). Suppose first that we replace r and 9 with r' and O' everywhere in 
the metric of equation (21). Clearly such a relabelling does not affect the solution. We 
then re-express the line element back in the original (r, 9) coordinate system. The result 
is that all of the fields (ai... d 2 ) have transformed to new fields. Four of these, 01 , a 2 , d\, 
d 2 transform according to the simple rule 


o-i(r , 9) a[(r, 0) = ai(r',9'). 


(30) 


Any field transforming in this manner is said to transform as a covariant scalar under 
diffeomorphisms. All of the 10 scalar functions introduced in the uii of equation ( |27| ) 
behave as covariant scalars. 
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The remaining variables in the g ^ vectors, b\, b 2 , ci, C2, pick up derivatives of the 
coordinate transformation under the diffeomorphism described above. The transformed 
variables are given by 

fb[(r,0) d 2 (r,6) \ _ fb\{r',0') c 2 (r',9') \ (d r ,r d r >6\ 

V& 2 M) c' 1 {r,0)/r) \b 2 (r',9’) Cl (r', 0')/r') \d e >r d e ,Q) ' ^ 

It is simplest for our purposes to view the eight transformed variables a[ ... d' 2 as a new 
set of functions obtained from the original set by a gauge transformation. As such, the 
two sets are physically indistinguishable, and this gauge freedom must be fixed before we 
can write down a unique solution. The standard, metric-based approach to solving the 
Einstein equations usually attempts to fix this gauge freedom at the outset by restricting 
the terms in the metric. For example we could write the line element in the form j|] 

ds 2 = e 2a dt 2 - e 2 ^(# - u dtf - e 2 ^ dr 2 - e 2 ^ 3 d9 2 . (32) 


This leaves us with a set of 5 variables, with the freedom to fix the relation between fi 2 
and //3 by a further coordinate transformation. Our strategy (in keeping with that of the 
NP formalism) is to leave this gauge unfixed until much later in the calculation. This is 
possible because the coframe derivatives transform as 

L,"L' r = 6i(/,«') A + c 2 (r',9') A (33) 

with a similar result holding for Lg. Provided we formulate all equations in terms of 
coframe derivatives of covariant scalars, the entire structure will transform covariantly 
under diffeomorphisms. That is, the gauge freedom can be ignored until later in the 
problem. We then find that certain physical fields emerge and it is sensible to equate these 
with combinations of the chosen coordinates. This is the point at which the diffeomorphism 
gauge freedom is fixed. Working in this manner ensures that the chosen coordinates have 
a direct physical interpretation. 

The second gauge freedom is that of applying a local Lorentz rotation to the g^ vectors. 
This is the gravitational analogy of a phase transformation in quantum theory. Under a 
Lorentz rotation the g^ vectors transform according to 

g 11 ~ Rg^R (34) 

where R is a local rotor. The reciprocal vectors g ^ transform the same way, and the metric 
derived by equation ( |20|) is therefore invariant. Clearly this freedom is a gauge symmetry, 
as it does not change any physical quantity. This symmetry is removed if one works directly 
with terms in the metric, but our approach is to keep the symmetry explicit and use it to 
simplify our equations. The u>i bivectors represent the connection for the gauge symmetry 
of equation (|34|). The combination 

Dti = + g fJ ,- r y l Ui>< (35) 

defines the covariant derivative for objects transforming in the manner of equation (|34j). 
For our axisymmetric setup the degrees of freedom in the Lorentz gauge are reduced from 


six to two. These are a rotation in the la^ plane and a boost in the a^ direction. The 
rotors describing these can be combined to form the single rotor 

R = ex.p(aa t p/2) exp(/3Ia t p/2) = exp(wIa t j ) /2 ), (36) 

where w = (5 — la is an arbitrary function of (r, 6). The freedom to apply this rotor will 
be used to simplify the Riemann tensor. 

One final (non-gauge) freedom present is that the equations are unchanged if all vari¬ 
ables are scaled by a constant amount, a± ca\, G cG, etc. This freedom is often 
employed to fix the asymptotic behaviour so that the vectors map to a flat space coor¬ 
dinate frame at spatial infinity. 

4 The Riemann tensor and vacuum fields 

The essential covariant object to construct is the Riemann tensor. The information con¬ 
tained in this is compactly summarised in a set of six bivectors. These are defined by 

Ry — L^LUj Ljijji —|— (jjj x tijj C'ij (jjfe. (37) 

Calculating each of the terms Ry is simply a matter of organisation, and is well suited to a 
symbolic algebra package. A useful set of intermediate results is provided by the identities 

L t a r + u> t x a r = (S + IK^Ia^ L t a g + ay x oy = —(T + IJ)Ia y 

L r a r + LU r x ay = (S + IK)a g L r a g + t o r x a g = —(S + IK)a r 

L g a r T loq xay = {G T IJ^)a g Goa q T ujq x a g = — (G 1J^a r 

L'par + u^xar = (G + IJ)a,p L^ao + u^x a e = (H + liTja^. (38) 
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On substituting the form for the 07 into the definition of Ry we now obtain the Riemann 
tensor in terms of abstract derivatives as 

R rt = [~L r (T + IJ) + (S + IK)(S + IK) + T(T + IJ) + I(I< + K)(H + IK)^ja r 

+ (~L r (S + IK) - (S + IK)(T + IJ) + T(S + IK) - I{K + K)(G + I J)) a 0 

R 0t = (~L g (S + IK) — (G + U)(T + IJ) + S{S + IK) + I(J - J)(G + IJ))a e 

+ ( -L e (T + IJ) + {G + IJ)(S + IK) + S{T + IJ) - I{J - J)(H + IK)\a r 

R # = -((G + IJ)(T + IJ) + (S + IK)(H + IK)^ 

RrO = - ( L r {G + 13) — Lq(S + IK) + G(G + 13) + S(S + IK )) la$ 

R r </, = [L r {G + IJ) - (S + IK)(H + IK) + G(G + IJ) + I(K - I<)(S + IK)\ lag 

- ( L r (H + IK) + {S + IK)(G + /J) + + IK) - I(K - K)(T + /J)) Ia r 

R 00 = ~(Lo(H + IK) + (G + IJ){G + IJ) + H{H + IK) + I(J + J)(T + /J))/u r 

+ (L e (G + IJ) - (G + IJ){H + IK) + H(G + IJ) - I(J + J){S + IK^Iag. 

(39) 

The standard coordinate expression of the Riemann tensor is obtained from the preceding 
bivectors by forming the contraction 

R^upa = (g^hg^-Rijj'-gai^gp- (40) 


As with the Christoffel connection, these coefficients cannot be written down until an 
explicit form is obtained for the g^ vectors. 

By working with a general form for the spin connection we have no guarantee that the 
Riemann tensor satisfies all of its required symmetries. The first of these to consider arises 
from the assumption that there is no torsion present. In this case the Riemann tensor 
satisfies 

1 i AR ij = l(fR ij + R ijl i ) = 0. (41) 


This says that a trivector vanished for each value of the index j, so corresponds to a set of 
16 constraints. The Riemann tensor is encoded in six bivector-valued functions, giving a 
total of 36 scalar degrees of freedom. Equation (41) reduces this to the familiar 20 degrees 
of freedom of a general Riemann tensor. 

In this paper we are interested in vacuum solutions to the field equations. For these 
the contraction of the Riemann tensor must also vanish, 


7 i xRi i = K7 < Rii-Rij7 i )=0. 


This can be combined with equation 


(fll|) to obtain 


7 ? Rij = 0 . 


(42) 


(43) 
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This equation is satisfied by all Weyl tensors. Equation ([£|) reduces the number of degrees 
of freedom in the vacuum Riemann tensor from 20 to 10 — the expected number for vacuum 
solutions. The compact combination of the symmetry and contraction information in 
equation |43| is unique to the Clifford algebra formulation. 

Before proceeding, it is useful to adopt a slightly different notation for the Riemann 
tensor. We define a linear map from bivectors to bivectors by writing 

R(B) = iB.(VA 7 ! )R li . (44) 

Equation ( |39| ) then shows that the Riemann tensor has the general form 

R(o>) = aqa> + (3\Oq R(/o>) = a 4 Ia r + p 4 Iag 

R(o- e ) = a 2 a e + /3 2 o> R (I erg) = a 5 Iag + fi 5 Ia r 

R(o>) = ct 3 cR (io>) = (45) 


where each of the a* and /3i are scalar + pseudoscalar combinations. We can now clearly 
see how a complex structure emerges based on the spacetime pseudoscalar. On writing 
out of the four equations (43) in full, and pre-multiplying each equation by 7 j, we obtain 
the four equations 


oyR(oy) + (JgR{ag) + o>R (o>) = 0 
ayR(oy) - IagR(Iag) - Icr ( t > R(Ia ( p) = 0 
—/oyR(Iay) + crgR(ag) - Ia^Rilcr^) = 0 

—Ia r R(Ia r ) - Ia e R(Ia e ) + a^a^) = 0. (46) 

Summing the final three equations and employing the first we obtain 

Ia r R(Ia r ) + Ia e R(I(j e ) + Ia^Ia^) = 0. (47) 

Substituting this back into the final three equations we obtain the simple relation 

R(M) = IR(vi), (* = r , 9, (j)). (48) 

This equation shows that the vacuum Riemann tensor (and the Weyl tensor in general) is 
linear on the pseudoscalar. This is the origin of the complex structure often employed in 
analysing the Weyl tensor. Equation (^) then sets 


= a 4 , a 2 = a 5 , a 3 = a 6 , = (3 4 (3 2 = f3 5 . (49) 

Equation (^7|) says that, viewed as a complex linear function, the vacuum Riemann tensor is 
symmetric and traceless. The most general form of R(-B) allowed for vacuum axisymmetric 
solutions therefore has the form 


R(o>) = aqoy + (3a g 
R(o-0) = a 2 a g + (3a r 

R(a>) = — (aq + a 2 )a ( j > . (50) 
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We can simplify equation (50) further by recalling the Lorentz gauge freedom present 
in our setup. Under a Lorentz gauge transformation parameterised by the rotor R the 
Riemann tensor transforms according to 


R (B) ^ R '(B) = RR(RBR)R. 


(51) 


In our axisymmetric setup the rotor R. is restricted to the form of equation (36). Under 
this transformation we find that (3 transforms to 


/3' = cos(2 w)/3 — ^(ai — a 2 ) sin(2u;) 


(52) 


where the scalar + pseudoscalar combination w is treated as a complex number. We can 
therefore eliminate the off-diagonal term from R(B) by setting tan(rc) = 2/3/(a± — 0 : 2 ). 
This fixes w and so removes the rotational degrees of freedom. The conclusion is that the 
Riemann tensor for the vacuum outside an axisymmetric source can be written 


R(oy) = a\<j r 
R(cr 6l ) = a 2 a e 

R (o>) = -(«i + a 2 )(74>- (53) 


This can be compared with equation (39) to obtain a set of 20 equations, which naturally 
couple into 10 scalar + pseudoscalar combinations. Four of these equations return explicit 
formulae for the action of L r and Lg on G , T, J, S, H, K, 


L r (S + IK) = -(S + IK) (T + IJ) + T(S + IK) - I(K + K)(G + IJ) 

L r (H + IK) = -(S + IK) (G + IJ) - G{H + IK) + I(K - K)(T + IJ) 

L e (G + IJ) = (G + IJ)(H + IK) - H(G + IJ) + I(J + j)(S + IK) 

L e (T + IJ) = (G + IJ)(S + IK) + S(T + IJ) - I(J -J)(H + IK). (54) 

A further four equations return the remaining coframe derivatives of G , T, J, S, H, K, 
but also include terms from the Riemann tensor 

L r (G + IJ) = (S + IK)(H + IK) - G(G + IJ) - I(K - K)(S + IK) + a 2 
L r (T + IJ) = {S + IK)(S + IK) + T(T + IJ) + I(K + K)(H + IK) - ai 
L e {S + IK) = -(G + IJ)(T + IJ) + S(S + IK) + I(J - J)(G + IJ) - a 2 
Lg(H + IK) = —{G + IJ)(G + IJ) - H(H + IK) - /(J + J)(T + IJ) + ai . (55) 

The final equations give 

ai + a 2 = (G + IJ) (T + IJ) + (5 + IK)(H + IK), (56) 

and 

L r (G + 13) -L g {S + IK) + G(G + 13) + S{S + IK) = ai + a 2 . (57) 

So far we are some way short of a fully determined system. The next step is to impose the 
Bianchi identities. 


12 





5 The bracket structure and Bianchi identities 


The Lie bracket structure of equation ( |28|) gives rise to a series of higher order constraints. 
This information is summarised in the Bianchi identities. For vacuum solutions the Bianchi 


identities take the simple form [10] 


7* Li(R(B)) - R (LiB) +u> t x R(B) - R(^ X B) = 0, 


(58) 


which holds for any bivector B. The linearity on / and the properties of the Riemann 
tensor imply that we only obtain new information for B = a r and a At this point it is 
useful to set 


a = a.\ + «2 

<5 = ai + 2a 2 - (59) 

The reason for this choice is that <5 = 0 corresponds to a type D solution. Applying 
equation ( |58|) to B = 07 , yields the pair of equations 

L r a = -(3 a - S)(G + IJ) + <5(T + IJ) 

Lga = (3o - 5){S + IK) - 5(H + IK). (60) 

These are entirely consistent with equations (|54|), ( [55| ) and (^) so contain no new infor¬ 
mation. Applied to B = a r , however, the Bianchi identity does provide two new equations, 

L r 5 = -(3ct - <5)(G + IJ) + 25(T + IJ) + (3a - 2<5)(G + 13) 

Lg5 = (3a -5)(S + IK) - 25{H + IK) - (3a - 2<5)(5 + IK). (61) 

We now have expressions for many of the coframe derivatives of the main physical 
variables. These derivatives must all be consistent with the bracket structure, which 
reduces to the single identity 


[L r ,Le\ = -SL r -GL e . (62) 

This identity already holds for all of the unbarred fields. The identity can also be applied 
to 5 and yields one further equation, 

(3a - 25) (L r (S + IK) + L e {G + IJ) + G(S + IK) - S{G + /J)) = 

3a((S + H + 2IK)(G + IJ - G - IJ) - (G + T + 2 IJ)(S + IK - S - IK )). (63) 

But this is as far as the equations can be developed without further physical information. 
There are no further expressions on which to evaluate the bracket identity, and we do not 
have a complete set of coframe derivatives. That is, the vacuum structure is currently 
under-determined. This is to be expected as there is no unique vacuum solution outside 
a rotating star. To obtain a unique solution we must either impose suitable boundary 
conditions, or make a further restriction on the form of the Riemann tensor. For example, 
we could set J, J, K and K to zero. The effect of this is to set a and 5 to real variables, so 
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that the equations reduce dramatically. The equations then describe the fields outside a 
static, axisynunetric source. That is, a non-rotating lump of matter which axial symmetry. 
This setup was first discussed by Weyl in 1917 [|l5[ ||. But before we embark on solving 
the general equations for a particular system, we give a comparison of our equations with 
the NP formalism. 


6 The Newman—Penrose formalism and complex structures 


In the present scheme the field equations for axisymmetric vacuum fields are summarised 

together with the bracket identity 


by equations (p4|) , (p5|), (p6|), (|57|) (|6lD and 
There is clearly a close analogy between these equations and the NP formalism. The first 
point to note is that all equations now consist of scalar + pseudoscalar combinations. As 
all other algebraic elements have been removed, the only remaining effect of the Clifford 
pseudoscalar I is to provide a complex structure through the result (§J) that 1 2 = —1. We 
can therefore systematically replace / by the unit imaginary i , 


/ i—> i. 


(64) 


This is quite helpful typographically, as it distinguishes field variables from the the (con¬ 
stant) pseudoscalar. We therefore employ this device at various points. It must be remem¬ 
bered, however, that when we construct the full solution in terms of the Wj, the imaginary 
i must be replaced by the pseudoscalar. 

Clearly, part of the complex structure in the NP formalism has its geometric origin 
in the properties of the spacetime pseudoscalar. This unites the complex structure with 
spacetime duality. But the starting point for the NP formalism is a complex null tetrad. 
We define the four (Minkowski) null vectors 


1 = Tf (7t + 7 ^ 

n = - 7 *) 


m = ^-={ rfr + i l0 ) 
fh = ~^= (q'r - ijo). 


(65) 


If these vectors replace the { 7 ^} in equation (19) the result is a complex null tetrad. The 
differential operators defined by the null tetrad are, following the conventions of Kramer 
et al. [Hi, 


1 


D — -^=(£t + A/>) 
A = -/=(Tt - L^) 


5 — —j=(L r + iLg ) 
5* = -L (L r -iLg). 


( 66 ) 


One can form bracket identities on these derivatives to obtain each of the spin coefficients 
in terms of the fields in the 07 . For example, we find that 


[A, D] = (7 + 7 *)D + (e + e*)A - (r* + tt)< 5 - (r + 7 r*)<T = [L t , L#] = 0, (67) 
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where 7 , e, r and n are spin coefficients. Continuing in this manner we obtain expressions 
for each of the spin coefficients. These include the relations 


r = 


2y/2 


( T-G + i(S-H )) 


and 


v = — 


2\/2 


(G + T + 2A' + i(2J- 5- A)). 


( 68 ) 


(69) 


There are two problems with the NP formalism. The first is that the key variables (in the 
spin coefficients) are inappropriate combinations of the natural variables in the 07 . The 
second is that there are two distinct complex structures at work. One is geometric, and 
has its origin in the pseudoscalar; the other is purely formal and is inserted by hand in the 
null tetrad. If both of these are represented by the same unit imaginary i then the true 
geometric structure is lost. This is a fundamental flaw, which compromises our ability to 
solve the equations. 

As an aside, it is possible to formulate the null tetrad ( |65| ) in terms of a real Clifford 
algebra. In this case the entire complex structure is provided by the pseudoscalar. A 
consequence of this is that the null tetrad consists of vector and trivector combinations (H| . 
This is appropriate for some encodings of supersymmetry, but is not suitable for analysing 
the gravitational equations, where the null tetrad must consist of complex vectors. 


7 Type D vacuum solutions 

In order to take our method on to a solution we need to impose some extra constraints 
on the vacuum fields. Ideally one would like to impose boundary conditions on some fixed 
2 -surface representing the edge of a rotating source, and then propagate the fields out into 
the exterior region. Here we follow a slightly different approach and restrict the algebraic 
form of the Riemann tensor to type D. The relationship between this restriction and the 
existence of a horizon is described in section (|l3f). 

Viewed as a complex linear function the eigenvalues of the Riemann tensor (^3|) are ctq, 
ctq and — (aq + aq). This reduces to a type D tensor if two of the eigenvalues are equal. One 
way of achieving this is to set ctq = 0 . 2 - But for this case the equation structure collapses 
to a simple system which is not asymptotically flat. The remaining two possibilities are 
gauge equivalent, so we can restrict to type D by setting 


Ctq = — Oq — Ctq. 

It follows immediately that <5 = 0. The Riemann tensor now reduces to 

R(o>) = 2acr r 
R(a e ) = -aoQ 
R(o>) = 


(70) 


where 


a = (G + IJ)(T + IJ) + (5 + IK)(H + IK). 


(71) 

(72) 
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The Riemann tensor can be encoded neatly in the single expression 


R (B) = ^ a(B + 3a r Bcr r ), 


an expression which is unique to the Clifford algebra formulation. 

Returning to equation (61) and setting 6 = 0 we see that we must now have 


G + IJ = G + IJ 
S + IK = S + IK. 


( 73 ) 


(74) 


This simplification for type D fields explains our notation using barred and unbarred 
variables. With this simplification we find that the derivative relations (54) and (|55| ) 
collapse to give 


L r (G + iJ) = -(G + iJ) 2 - T(G + iJ) 

L r (T + iJ) = (S + iKf - (2 [G + iJ) - T) (T + iJ) - 2 S(H + iK) 

L r (S + iK) = —iJ(S + iK) - 2 iK{G + iJ) 

L r (H + iK) = —{G + iJ)(S + iK) - G(H + iK) (75) 


and 


L g (S + iK) = (S + iKf + H(S + iK) 

L e (H + iK) = —(G + iJ) 2 + (2(5 + iK) - H){H + iK) + 2 G(T + iJ) 

L e (G + iJ) = iK(G + iJ) + 2 iJ{S + iK) 

Lq(T + iJ) = {G + iJ)(S + iK) + S(T + iJ ). (76) 

We have here adopted the convention of representing the pseudoscalar I with the unit 
imaginary i. This is a useful device when analysing this type of equation structure with 
a symbolic algebra package. The above equations are all consistent with the bracket 
structure, which now reduces to the single identity 


[L r ,L e \ = -SL r -GLg. 


(77) 


This set of equations is now complete — we have explicit forms for the intrinsic derivatives 
of all of our variables, these are all consistent with the bracket structure, and the full 
Bianchi identities are all satisfied. Obtaining a set of ‘intrinsic’ equations such as these is 
a key step in our method. 

The equations (75) and (76) exhibit a remarkable symmetry. The L r and Lg derivatives 
can be obtained from each other through the interchange 


Gh5 J^K 

T <-> H L r — Lg. (78) 

The origin of this symmetry is explained in Section (^), where it is related to the conjugacy 
transformation discussed by Chandrasekhar Q. 
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8 The Schwarzschild solution 


Before proceeding to the Kerr solution, it is helpful to look at how the Schwarzschild 
solution solution fits into our scheme. The restriction to spherical symmetry sets all of 02 , 
62 , C 2 and c ?2 to zero, and we also require that d\ = c\. The remaining functions, a±, b\ and 
ci are all functions of r only. The bracket relations now tell us that the only remaining 
functions in are T, G = G and H, with H given by 


ci cos 9 
r sin# 

We are therefore left with the pair of equations 


(79) 


L r G = -G 2 - GT 

L r T = -2 GT + T 2 , (80) 


where L r = b\d r . Dividing these, and solving the resulting homogeneous differential 
equation, we obtain 

GT = A 0 (G 2 - 2GT) 3 / 2 , (81) 

where Ao is an arbitrary constant. A significant point is that this relation is derived without 
fixing the diffeomorphism gauge. The relationship is an intrinsic feature of the solution. 
The Riemann invariant a is now simply GT, and satisfies 


L r cx — —3 Got. 


(82) 


The equation for LgH reduces to the algebraic expression 


-T = G 2 - 2 GT. 

(83) 

8 

a = GT = A,; — . 

(84) 


If follows that 


We retain the diffeomorphism freedom to perform an r-dependent transformation, which 
we can employ to fix the functional dependence of one of our variables. The obvious 
way to employ this freedom now is to set ci = 1. This ensures that r defines the proper 
circumference of a sphere in the standard Euclidean manner. The gauge choice also ensures 
that the tidal force, controlled by a, falls off as r~ 3 , which agrees with the Newtonian result. 
It is then clear, from comparison with the Newtonian result, that Ao can be identified with 
the mass, so 

M 

a = --j. (85) 

This ensures that the status of r is lifted from being an arbitrary coordinate to a physical 
held controlling the magnitude of the tidal force. 


Equation (82) now tell us that G is given by 


( 86 ) 
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and equation (|^) becomes 


(87) 


This integrates to give 


1 dr \r ) r 2 r 3 ’ 


b\ 2 = cq — 2 M/r 


( 88 ) 


where Co is the arbitrary constant of integration. A combination of a further constant 
rescaling of r and a global scale change can be used to set co = 1, so that g r tends to a 
flat space polar coordinate vector at large r. 

Finally, we recover a\ from the \Lt,L r \ bracket relation, which yields 


L r a i — Tax- 


(89) 


Since L r b\ = — Tb \, it follows that ai&i must be constant. This constant can be gauged 
to 1 by a constant rescaling of the time coordinate t, so we have 


ai = (l-2M/r)' 1/2 . 


(90) 


This solves for all of the terms in {g^}, and calculating the associated metric recovers the 
Schwarzschild line element. 

This derivation illustrates a number of the key features which reappear in the derivation 
of the Kerr solution. The intrinsic differential structure can be taken almost completely 
through to solution without ever introducing any coordinatisation. Coordinates are only 
fixed in the final step, when one solves for the metric coefficients. For the Schwarzschild 
case this results in equation ([87]), which is solved directly be integration and could have 
been computed numerically had a simple analytic solution not been available. All of the 
remaining terms are found algebraically, or by integration alone. Another feature is the 
employment of any residual gauge freedom to set the asymptotic conditions so that the 
{d/i} vectors map to a flat space coordinate frame as r i—> oo. Finally, the form of the 


solution shows that the fields are only valid for r > 2 M. Our initial ansatz (17) is not 
sufficiently general to describe a global solution when a horizon is present. This problem 
is only resolved by introducing further terms to allow for a t-r coupling [1C]. 


9 Simplifying the Kerr equations 

We are now in a position to complete the derivation of the Kerr solution. Our approach is 
to construct a solution which is asymptotically flat at large r, for all angles. This condition 
turns out to be sufficient to construct a unique solution, which is then easily seen to be 
the Kerr solution. 

The first step in solving the equation structure of (ff^) and ( ff6| ) is the identification of 
a number of integrating factors. If a pair of variables A and B satisfy the equation 

L e A - L r B = GB + SA (91) 

then an integrating factor x exists defined (up to an arbitrary magnitude) by 

L r x = Ax , Lgx = Bx. (92) 
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Equation (^) ensures that the definition of the integrating factor is consistent with the 
bracket structure ([f7|). This means that, when a coordinatisation is chosen, the equations 
for the partial derivatives of x are consistent, and x can then be computed by straightfor¬ 
ward integration. 

The first integrating factor is constructed from a. The factor of 3 in equation ( |60| ) 
prompts us to define 

Z = Z 0 a~ l/3 , (93) 

where Zq is an arbitrary complex constant. From equation (|60|), Z satisfies 

L r Z = (G + iJ)Z 

LgZ = -{S + iK)Z. (94) 

We can therefore use Z to simplify a number of our equations. Separating into modulus 
X and argument Y, 

Z = Xe iY i (95) 

we find that 

L r X = GX 

L g X = -SX, (96) 

which provides an integrating factor for G and S. This is particularly important, as a 
comparison with the bracket relation (|T7|) now tells us that 


[XL r ,XL g \ = 0. (97) 

Recovering a pair of commuting derivatives like this tells us that we should fix our dis¬ 
placement gauge freedom by setting 62 = C 2 = 0. With this done, we can write 


XL r = b(r)d r 

X Lg = c(0)dg, (98) 


where b(r) and c(9) are arbitrary functions which we can choose with further gauge fixing. 

A search through the remaining derivatives reveals that the pair T and — H satisfy an 
equation of the form of (91). We can therefore introduce a further integrating factor, F, 
with the properties 

L r F = TF, LgF = -HF. (99) 


Both Z and F are unchanged by the conjugacy operation of equation (78). 

Now that we have suitable integrating factors at our disposal, we can considerably 
simplify our equations for G + i J and S + iK to read 


L r (FZ{G + iJ )) = 0 
L e (FZ(S + iK )) = 0 
L r (XZ(S + iK)) = 2 XZ{SG + JK) 

L e ( XZ[G + iJ)) = -2 XZ{SG + JK) (100) 
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These equations retain the symmetry described by equation (78). Equations ( |10C| ) now 
focus attention on the quantity SG + JK. Computing the derivatives of this, we find that 


L r (SG + JK) = —{G + T)(SG + JK) 
L e (SG + JK) = (H + S)(SG + JK). 

It follows that 

L r (XF(SG + JK)) = L e (XF{SG + JK)) = 0, 


( 101 ) 


( 102 ) 


and hence that XF(SG + JK) is a constant. It turns out that the equations can be solved 
if this constant is non-zero, but the solutions require infinite disks as sources and are not 
asymptotically flat in all directions. These solutions are presented elsewhere. Since we are 
interested here in asymptotically flat vacuum solutions we must set this constant to zero, 
which implies that 

SG + JK = 0. (103) 


Equations ( 100|) now yield 


L r (XFZ 2 (G + iJ){S + iK)) = L e (XFZ 2 {G + iJ)(S + iK)) =0. (104) 


It follows that 


XFZ 2 (G + iJ)(S + iK) = Cl, 


(105) 


where C\ is an arbitrary complex constant. 

Remarkably, we are now in a position to separate the remaining equations into radial 
and angular derivatives. With L r and Lg defined by equation (98), we see that we can set 



FZ(G + iJ) = W(9) 

FZ(S + iK) = U (r). 

(106) 

We therefore have 


XW{6)U{r) = C\F, 

(107) 

and 

XZ(G + iJ) = C\/U[r) 

XZ{S + iK) = C\W{6). 

(108) 


This provides a structure consistent with our derivatives of G + i J and S + iK. 
We next form 


W(0) _ G + iJ _ .SJ - GK 
U(r) ~ S + iK ~ 1 S 2 + K 2 ’ 


(109) 


which is a pure imaginary quantity. This implies that W and U are n/2 out of phase. 
Since these are separately functions of r and 9, their phases must be constant. Returning 
to the derivatives of Z we see that 


XL r Z = b(r)d r Z = XZ(G + iJ) = Ci/R(r) (110) 

and 

XL e Z = c(0)d e Z = XZ(S + iK) = Ci/W(0). (Ill) 
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It follows that Z must be the sum of a function of r and a function of 0. Furthermore, 
these functions must also have constant phases ir/2 apart. Since the overall phase of Z is 
arbitrary (Z was defined up to an arbitrary complex scale factor), we can write 

Z = R(r) + M(0) (112) 

where R(r ) and ’!'(#) are real functions. With Z reduced to this simple form we can solve 
the remaining equations by finite power series. 


10 Power series solution 


With Z given by equation m , we now have 

XZ(G + iJ) = XL r Z = XL r R = R' 


and 

XZ{S + iK) = —XLqZ = -iXL e m = -i 4’. 
Here we have introduced the abbreviations 

R! = XL r R 
T = XL e ^>. 


(113) 

(114) 


(115) 


These help to encode the fact that XL r R is a function of r only, and XLq^> is a function 
of 6 only. The functions implied in the R' and notation are not necessarily pure partial 
derivatives, though they can be chosen as such with suitable gauge choices. 

We now have 

F _ i 
X ~ CiR’V' 

which implies that C\ is an imaginary constant. Since F is only defined up to an arbitrary 
scaling, we can choose this constant as i and write 



F _ 1 

x ~ ~m>' 


Differentiating this relation now yields 


X(T-G) = X^L r ^ = 


XLrR' 

R' 


RT_ 

~w 


(117) 


(118) 


and 


X(H - S) 


X F XLg'F T 
X F Le X ~ T _ T 


(119) 


where we have introduced the obvious notation R" and T for the second derivatives. 

The final relation we need to satisfy is to construct a from R and and equate this 
with ( Zo/Z ) 3 . On writing 


a = (G + iJ) 2 + (T - G)(G + iJ) + {S + iK) 2 + (H - S){S + iK) (120) 
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we see that 


X 2 Z 2 a = Ft! 2 - R"(R + M) - T 2 - i<S{R + ) 


= R! 2 - RR” - T 2 + mi - i(R+ VR). 


4,2 


This must be equated with 


X 2 Z 2 a = Z 0 6 X 2 /Z = ZtfZ* = Z q 6 [R - iT) 


( 121 ) 


( 122 ) 


with Zo 3 an arbitrary complex constant. These are the only differential relations we need 
to satisfy in order to obtain a valid solution. 

We could proceed now by fixing the remaining displacement gauge freedom and speci¬ 
fying forms for R and \k, or d{r) and c(9). While this approach certainly works, we can in 
fact do better by continuing to work in our abstract, intrinsic fashion. The key is to notice 
that R' 2 and T 2 must be power series in R and T respectively in order for the right-hand 


side of equation (121) to be linear in R and \k. We therefore set 


R' 2 = k2 R 2 + k\R + ko 
T 2 = / 2 ^ 2 + ^id' + / 0 , 


(123) 


where the ki and l t are a set of six constants (any higher order terms vanish). These 
equations are differentiated implicitly to yield 


R" = k 2 R + k 1 /2 

T = l 2 ^ + h/2. 


(124) 


We now substitute these series into equation (121). In order to remove the constant and 
terms we must set 

h = k 0 , h = -k 2 . (125) 


The remaining terms factorise to yield 


X' 2 Z 2 a = \{k x - ih)(R - M), 


(126) 


which now satisfies equation (122). We have therefore solved all of the equations ‘intrinsi¬ 
cally’. This process has revealed the presence of an arbitrary complex constant controlling 
the Weyl tensor, which must be understood physically. The remaining tasks are to pick 
a suitable coordinatisation and to integrate the bracket relations to reconstruct the g^ 
vectors and the metric. 


11 Integrating the bracket structure 

In arriving at an intrinsically-defined power series solution we have completed much of 
the work required to generate a solution. In particular, we can now integrate the bracket 
relations (28) to solve directly for the coframe vectors, and hence for the metric. In so 
doing we have to include a number of constants of integration. Some of these are gauge 
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artifacts and others can be removed by imposing either elementary or asymptotic flatness. 
The end result of this process, as we show below, is that the constant l\ must be set 
to zero to achieve a solution which is well-defined everywhere outside the horizon, and 
asymptotically flat. A choice of position gauge can then be employed to set 

M 

a (r — iL cos#) 3 


(127) 


This could then be used to find the coframe vectors directly, but it is more instructive to 
see how the bracket structure integrates intrinsically. 

The bracket relations (28) produce a series of relations for a\, a 2 , b\ and b 2 . By applying 
the known integrating factors these become 


L 

L r (x 


1 \ 2 K a 2 

02 Wo 


r ' F 1 / F r sin# 


L r IT 

hr 


r sin# / 
di 

r sin# 

1 


= 0 

2 K di 


Le 


L(){Xa\) — 0 
1 a 2 \ 2 J 


F 2 j F r sin# 


F r sin# J F 01 
1 di \ _ 2J 
F r sin# J F 2 

Lg(Xd 2 ) = 0. 


To simplify these equations (128) we first form 


and 


2 J 

- 

2 T'k 

XF 

X 2 

x 5 J 

2 K 

= = 
X 2 

2 RR' 

XF 

X 5 


R' = L e 


R' 2 \ 
X 2 / 


' 9 ( ® 


The set of equations (128) therefore integrate simply to yield 

a\R' a 2 4> 


X 

hi 


d 2 R' 


X 


Xr sin# 
d ,\T 
Xr sin# 


= <Si 


= ei 


(128) 


(129) 

(130) 


(131) 


where <5i and e\ are two arbitrary constants. But we know that Xaq and Xd 2 are inde¬ 
pendent of 9, and Xdi/(rsin#) and Xa 2 /(r sin#) are independent of r. We can therefore 
separate the above relations and write 


a 1 
d 2 


<5o T 5i R 2 
XR' 

eo + ei R" 
XR' 


a 2 _ So — di^ 2 
r sin# X T 

di _ e 0 - eiT 2 
r sin# X T 


( 132 ) 


where <5 q and £q are two further constants. 
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The next problem to address is the behaviour of the solution on the 2 -axis. The 
problem here is the ambiguity in the coordinate system at the two poles. The situation 
can be easily sorted out by converting back to Cartesian coordinates and demanding that 
the equivalent vectors {g *, g x , g y , g z } are well-defined on the axis. We quickly see that on 
the axis we must have <22 = 0 and d\ = c\. Similar considerations hold for the cu* bivectors. 
For these we see that T must vanish on the axis and so must contain a factor of sin#. It 
follows that 

l 2 ^(0) 2 + Zi'F(O) + l 0 = + l 0 = 0. (133) 

A further requirement from the u>i bivectors is that we must have 

eo — eiT 2 — cos#*F = 0, # = 0,7r (134) 

which tells us that T, and hence T, must change sign between 0 and ir. Since \F 2 must 
have the same value at 0 and n to ensure that <22 vanishes on the axis, we see that we must 
have 

T(0) = -’F(tt). (135) 

We cannot have T = 0 at the poles, so equations (|133| ) and ( |135| ) enforce l\ = 0. Solutions 
with a non-zero value of can be constructed, but these are only valid in the upper or 
lower half planes separately. Matching these together introduces an infinite disk of matter 
as a source, and the solutions are not true vacuum. These will be discussed elsewhere. 
With li =0, we now have 

'F = -k 2 ^. (136) 

A sensible gauge choice, which produces a well-defined solution for all 9, is to set c(#) = 1 
in equation (|9^), and scale the solution so that k 2 = 1. It then follows that 'F is a linear 
combination of sin# and cos#. Since T vanishes at the poles, we must have 

V(6) = —L cos#, l 2 = -1, l 0 = L 2 . (137) 

The Riemann tensor is now controlled by the complex function k\(R— iLcos #) -3 /2, which 
determines the magnitude of the tidal forces experienced in different directions. A large 
distance from the source we expect this to tend to the Schwarzschild value of — M/r 3 , so 
a sensible choice for the radial gauge is 


R = r, k\ = —2 M. 


From equation (98) it follows that 


b(r) = XL r r = R', 


(138) 


(139) 


and from equation (123) we have 

R' 2 = b{rf = r 2 - 2 Mr + L 2 = A. 

Here we have introduced the standard symbol A for the radial function b{r) 2 . 


(140) 
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With the preceding gauge choices we now have 


G + iJ = 


A 1 / 2 


T — G = — 


p(r — iLcosO) 
r — M 


pA 1 / 2 


S + iK = 
H — S = 


—iL sin# 
p(r — iL cos#) 
cos # 
p sin# ’ 


where 


2 y2 2 i t 2 2/i 

p = X = r + L cos 0. 


The equation for T shows that a horizon exists at A = 0. We also now have 

A 1 / 2 


Similarly, we now have 


a i = 


pA 1 / 2 


b\ = 


(<y 0 + hr 2 


p 


r 

Cl = 

P 


d 2 = — (5o ~ d'l L 2 cos 2 0) 
Lp 


d\ = — (eo — eiL 2 cos 2 0) 
Lp 


d 2 = 


1 


pA 1 / 2 


(e 0 + eir 2 ). 


(141) 

(142) 


(143) 


(144) 


In order that the g^ vectors approach a flat space coordinate frame at infinity we set = 1 
and e\ = 0. The former corresponds to a constant rescaling in r and the latter is achieved 
with a constant rotation, both of which are gauge transformations. The on-axis conditions 
d\ = ci and a 2 = 0 set eo = L and do = L 2 respectively. We then finally arrive at the 
vectors 


9 t = 


if = 


r 2 + L 2 t L sin0 


pA 1 / 2 

A 1 / 2 

-7 


■7* + 


</ = f" 

P 


9 r = 


-Y + 


L 


psind pA 1 / 2 


f- 


(145) 


The reciprocal set is given by 


A 1 / 2 L sin# 

9t = - 7 1 - 70 

P P 

9r = f±fj2 lr 

9e = Pie 

(r 2 + L 2 )sm9 LA 1 / 2 sin 2 # 

#0 = 2 - L -70- It- 146 

P P 

The metric derived from these vectors gives the Kerr solution in Boyer-Lindquist coordi¬ 
nates. 
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12 Conjugate solutions 


At various points in the text we have pointed out the symmetry in the equations described 
by equation (f7S|). We are now in a position to understand how this comes about, and 
relate the symmetry to the conjugation operation discussed by Chandrasekhar [[|]. We 
first need the general version of equation (128) applicable when the barred and unbarred 
variables are not necessarily equal. These are: 


di 


Lg 

Lg 


r r sin# 
«2 

r sin# 

di 

r sin# 
«2 


= -G 
= -G 
= -H 
= -H 


di 


r sin# 
Q2 

r sin# 
di 

r sin# 
02 


- (K - K)d 2 

- (K - K)ai 

+ (J 4- J)d 2 

+ («/ + J)ai 


Lga\ = Sai — (J — J) 


a 2 


r sin# 


Lgd 2 = Sd 2 - (J - J) dl 


L r a\ =Ta\ + (K + K) 
L r d 2 = Td 2 + (K + K) 


r sin0 
02 


r sin0 r sin@ 

The full conjugacy symmetry is therefore described by 

oi ^ —-7 d 2 

r smti 

L r —Lq T H 

Go S Go S 

J o-> K J ^ K. 


r sin@ 
di 

r sin0 


02 


r sin@ 


(147) 


(148) 


From the form of the Riemann tensor ([ll]) we can see that this interchange simply swaps 
round the various coefficients in R(-B). For a general, non-vacuum configuration this in¬ 
terchange will affect the Einstein tensor and therefore alter the matter distribution. For a 
vacuum solution, however, the Riemann tensor reduces to the form of equation (|53|). This 
is unchanged under the above transformation, which therefore does yield a new vacuum 
solution. 

In terms of the metric, we already know how the gt and g$ terms transform. The 
L r <-> —Lq interchange is achieved by 


b\ —b 2 /r, c 2 «-> — c\/r. (149) 

The effect of these interchanges on the line element (^) is to leave the dr 2 , dr d6 and dO 2 
coefficients unchanged, and to transform the remaining terms in the following manner: 

Adt 2 + 2 Bdtdcj) — Cdcj) 2 <-» Cdt 2 + 2 Bdtdc/) — Adcj) 2 . (150) 

This is precisely the conjugation operation described by Chandrasekhar j|]. This operation 
produces a new vacuum solution because the procedure can be viewed mathematically as 
the complex coordinate transformation t <-> i(f>. Our approach reveals more clearly the 
effect of the transformation on the various physical terms in the uji bivectors and the 
Riemann tensor. In particular, the interchange of G and S is likely to be problematic, 
as S is required to vanish on the 2 -axis, whereas no such restriction exists for G. This is 
certainly true of the Kerr solution, for which conjugacy symmetry does not yield a globally 
well-defined solution. 
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13 The Carter—Robinson Uniqueness Theorem 


One of the most remarkable features of the Kerr solution is that it represents the unique 
vacuum solution to the Einstein equations which is stationary, axisymmetric, asymptot¬ 
ically flat and everywhere nonsingular outside an event horizon. This result was finally 
proved by Robinson [jl7|, building on some initial work of Carter ]l 8 |. The proof proceeds 
in an analogous manner to standard uniqueness proofs by constructing an integral over 
the space outside the horizon. The boundary conditions are then employed carefully to 
show that a term representing the difference of two solutions must vanish [||. The proof is 
highly mathematical, so it is instructive to see how the present formalism highlights what 
is going on. 

So far, we have shown that the Kerr solution is the unique type D axisymmetric vacuum 
solution which is asymptotically flat and nonsingular outside the horizon. This proof is 
constructive — we constructed the most general solution with these properties and it 
turned out to be the Kerr solution. This is useful, as it shows that the role of the horizon 
is to force the Riemann tensor to be type D. To see how this comes about, we must return 
to the general set of equations before the restriction to a type D vacuum was made. These 
are embodied in equations (54), (|55l), 


and (63), together with the bracket 


identity (|62|). 

Our first task is to identify the horizon. The definition of the horizon is the surface 
over which the bivector defined by the Killing vectors gt and g^ is null. The equation for 
the horizon is therefore 

(gt/\g<p) 2 = o. (151) 

The volume element defined by the { g frame does not vanish anywhere, so an equivalent 
condition is 

(g r Ag 9 f = 0. (152) 

Since we are free to set C 2 = 62 = 0 with a choice of displacement gauge, the horizon in 
the present setup is characterised by g r2 = 0, or bi = 0. It follows that, at the horizon, L r 
acting on any finite, continuous quantity must vanish. 


From the \L r ,Lg\ bracket relation (28) we see that 


G = - L r 

ci 


(153) 


Since ci must be finite at the horizon, it follows that G vanishes there. We also have 


S= 1-Lob! = ^dgh. 

b\ rbi 


(154) 


But S cannot be singular at the horizon, as this would imply that the Riemann tensor 
itself was singular. It follows that b\ must take the form 


h = (r- r 0 ) v bi(r,6) 


(155) 


where ro defines the horizon, which is now a surface of constant r, 17 is some arbitrary 
exponent, and b± is finite at the horizon. 
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Next, we consider the force necessary to remain at a constant radius. This is simplest 
to analyse on the z axis. An observer remaining at rest on the z-axis must apply a force of 
magnitude T in the outward direction. This must be singular at the horizon, otherwise it 
would be possible to escape. This argument holds at all angles, and the horizon is therefore 
also defined (in this gauge) by T i—> oo. 

The invariants in the Riemann tensor, a and 5 must be finite at the horizon, so G and 
J must be zero there to ensure that the product T(G + iJ ) in a is finite. We also know 
that L r a and L r 8 must vanish at the horizon. Since T is singular there, equations ( | 6 l| ) 
force 8 = 0 at the horizon. It further follows from the equation for L r 8 that G + iJ must 
vanish at the horizon (and hence are equal to their unbarred counterparts). Finally, since 
the horizon is a surface of constant r, and 8 vanishes there, 8 must take the form 

8 = (r — ro) v 'S(r, 9), (156) 

where rj is a positive exponent and 8 is finite at the horizon. It follows that Lg8 must also 
vanish at the horizon, which in turn sets S + iK = S + iK there. 

The above argument shows that, at the horizon, the Riemann tensor is of type D and 
all the barred variables are equal to their unbarred counterparts. These are precisely the 
conditions met by the Kerr solution. In the light of this result, it is perhaps less surprising 
that the vacuum outside a horizon is described by the Kerr solution. Of course, this is not a 
complete proof of uniqueness, but it does serve to make the result more physically natural. 
One can take the above chain of reasoning further to show that the higher derivatives of 8 
also vanish at the horizon, from which it can be argued that 5 = 0 throughout the exterior 
vacuum. As shown earlier, this alone is sufficient to restrict to the Kerr solution. 

The main weakness in the preceding argument (which is shared with all proofs of the 
uniqueness of the Kerr solution) is the employment of a ‘bad’ gauge for the solution. The 
vectors and Ui bivectors are only defined outside the horizon, and various terms go 
singular on the horizon. To avoid this problem one must work with a more general ansatz 
for the gt 1 vectors which allows some coupling between t and r. Only then can one write 
down globally valid solutions. It would certainly be useful to have a version of the above 
argument which does not rely on a bad choice of gauge, and this is a subject for future 
research. 


14 Conclusions 


We have introduced a new framework for the study of axisymmetric gravitational fields. 
The technique draws on results from the gauge treatment of gravity, and builds on earlier 
work on spherical [1C] and cylindrical |Kj systems. In this paper we have concentrated 
on the vacuum equations. Further publications will detail how the technique is applied 
to matter configurations. We believe that this approach can help make progress in the 
long-standing problem of finding the fields inside and outside a rotating star. It is also 
expected that this framework should be helpful in analysing the disk solutions proposed by 
Pichon and Lynden-Bell p(l, and Neugebauer and Meinel [21]. The equations with matter 


included will be presented elsewhere. A significant feature of these is that a natural set of 
relations for the fluid velocity field emerge which can be carried into the vacuum sector [22], 
This velocity field in turn provides a physical basis for the Ernst equation. 
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The vacuum equations developed here are under-determined. They are only capable of 
unique solution when further physical restrictions are applied. The restriction to type D 
vacua is sufficient to yield a fully-determined set of consistent, ‘intrinsic’ equations. The 
development of such a set is similar to the process of obtaining an involutive set of equations 
in the formal theory of partial differential equations Q. The point of the latter process is 
that the equations can then be solved by a power series. In this respect, the power series 
solution of the equations developed here is remarkable. It hints at a power series solution 
method which does not rely on any definite coordinatisation. This idea warrants serious 
further investigation. 

We have shown that the twin restrictions of a type D vacuum and asymptotic flatness 
are sufficient to force us to the Kerr solution. The uniqueness theorem for the Kerr solution 
can therefore be interpreted as showing that the presence of a horizon forces a restriction on 
the algebraic form of the Riemann tensor. This provides a useful physical understanding 
of the uniqueness of the Kerr solution. Some indication of how this restriction comes 
about can be seen in the more general vacuum equations. But work remains to recast the 
formalism in a gauge which is not singular at the horizon. 

A future goal of this approach is to streamline and refine the techniques to the extent 
where time dependence can be included. It would then be possible to study the formation 
of rotating black holes, and their spin-up under accretion. Our approach opens up many 
new possibilities for the study of such systems. 
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